Roughness of moving elastic lines — crack and wetting fronts 
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We investigate propagating fronts in disordered media that belong to the universality class of 
wetting contact lines and planar tensile crack fronts. We derive from first principles their nonlinear 
equations of motion, using the generalized Griffith criterion for crack fronts and three standard 
mobility laws for contact lines. Then we study their roughness using the self-consistent expansion. 
When neglecting the irreversibility of fracture and wetting processes, we find a possible dynamic 
rough phase with a roughness exponent of = 1/2 and a dynamic exponent of z = 2. When including 
the irreversibility, we conclude that the front propagation can become history dependent, and thus 
we consider the value — 1/2 as a lower bound for the roughness exponent. Interestingly, for 
propagating contact line in wetting, where irreversibility is weaker than in fracture, the experimental 
results are close to 0.5, while for fracture the reported values of 0.55-0.65 are higher. 

PACS numbers: 62.20.Mk,68.08.Bc,05.40.-a 



I. INTRODUCTION 



The propagation of fronts or lines in disordered media occurs in a number of fields In particular, the propagation 
of cracks in heterogeneous materials and the wetting of disordered substrates involve lines with long-range elasticity. 
These two cases have attracted much attention over the past few years due to the difficulty to match theoretical 
predictions with experimental measurements. A commonly studied quantity is the so-called roughness exponent C, 
. which measures the roughness of the elastic line. For crack propagation, three different exponents must be distin- 
guished 0, one describing the roughness in the direction perpendicular to the crack propagation, a second the 
roughness in the direction of the propagation, and a third one (which will interest us in the following) describing the 
■ O in-plane roughness of the crack front during its propagation through the material. For the spreading/retraction of a 
C liquid wedge on a flat substrate, the contact line has only one roughness exponent. If h{x, t) stands for the transverse 
^ perturbation of the position of the line, the linearized equation of motion takes the non-dimensional form 

dh , . f h(x' ,t) — h(x,t) dx' 

> ', 

^ for both tensile cracks [1,0,0] and contact- lines 0,11, Here and elsewhere VV stands for the principal value of 
the integral. Solutions of stochastic growth models exhibiting scaling behavior are described by the time dependent 
N . correlation function 

^: {[h{x,t)-h{x\t')f)^\x-xr<f(^-^^^ , (2) 

!>■ '. 

■ where C is the roughness exponent of the interface and z is the dynamic exponent. From now on, the brackets (• ■ ■ ) 
^ [ denote average over disorder. However, for a medium with uncorrelated heterogeneities, the roughness exponent 
■'-j . corresponding to Eq. (jT]) was found d, 0, [l3l to be ^ = 0.39, whereas experiments yield ( = 0.55 to 0.65 for tensile 
/\ ' cracks and slightly more than C = 0.5 for contact lines [T^. Il5l|. 

To resolve this contradiction, besides accounting for additional physical effects (e.g. [IB])) it was proposed [itI [TsI. 
HBj [13) HH to include quadratic terms in h{x,t) in the equation of motion. However, different roughness exponents 
were predicted using the one-loop renormalisation group-^ = 0.45 [3, the self-consistent expansion-i^ =1/2 or 
numerical simulations-^ = 0.37 21]. When looking at the exact form of the equations used in these studies, it turns 
out that nonlinear terms have different forms. This is the first motivation of thepresent study, i.e. establishing the 
nonlinear equations of motion on firm grounds. In addition, experiments [HI) [13) O) 0) [iBl are mostly interpreted [^, 
IB) E) [3 in terms of the dcpinning transition of the front. Underlying assumptions are that the driving force is 
approximately constant, and that the statistical properties of the shape of the front are imposed at its dcpinning 
threshold. In fact, the front can achieve high velocities in experiments, well above the dcpinning threshold. Besides, 
the recent numerical simulations |,21|] suggest that the roughness at dcpinning is not affected by nonlinear terms in 
the equation of motion. Therefore, we will investigate the roughness of a front advancing at finite velocity. 

Here we study the dynamic roughness of lines advancing according to the nonlinear equations that we obtained. 
We find a rough phase characterized by an exponent C = 1/2, close to that observed for contact lines in wetting and 
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slightly below those measured for crack fronts. We suggest that the experimental roughness results from the freezing 
of an advancing rough front, so that it should be not only close to C = 1/2 predicted here, but it should also be history 
dependent, a fact that could be checked in experiments. We will use the Self-Consistent Expansion (SCE), a method 
developed by Schwartz and Edwards [13, and applied successfully to the Kardar Parisi-Zhang (KPZ) equation ^23|. 
The method gained much credit by being able to give reasonable predictions for the KPZ scaling exponents in the 
strong-coupling phase above one dimension where many renormalization group (RG) approaches failed [25l |. Another 
point which is especially relevant for our purpose is that for a family of models with long-range interactions of the 
kind treated here, SCE reproduced exact one-dimensional results while RG failed to do so (26j . 

The article is organized as follows. In section II, we derive from first principles the nonlinear equations of motion 
for crack fronts, and for the three available models for contact lines. We set all equations into a single framework and 
then introduce the annealed noise approximation corresponding to a line advancing at finite velocity. In section III, 
we describe the Sclf-Consistcnt Expansion and apply it to the nonlinear equation of motion, in order to investigate the 
possible phases and their roughness. Finally, in section IV we compare our results with previous works and discuss 
the relevance of our approach to experiments. 



II. EQUATION OF MOTION OF AN ELASTIC LINE 



A. Fracture 



1. The stress-intensity factor 




FIG. 1: Schematic representation of a tensile crack in an infinite body propagating in the plane y = Q. The average penetration 
of the crack front in the x-direction is L. The straight reference front in the z-direction and the perturbation h{z, t) around it 
are also shown. 



In order to write the equation of motion of a crack front, we first recall the form of its stress intensity factor P.[l9|. 
Consider a crack submitted to a tensile loading and propagating in the (x, z) plane, with a front located at a; = t), 
t being the time coordinate (see Fig.[T]). The static stress intensity factor K[h] was computed perturbatively to second 
order in h and is given by [l9[ 



which might be rewritten as 



\h'^{z)+VV 



h{z') - h{z) 
{z' -zf 



l + VV 



hjz") - hjz') dz" 
{z" - z'Y 



Az!_ 
27r ' 



Kq ' ' 8 4 



n[hn[h]] + o{h'^) 



(3) 



(4) 



Here is the stress intensity factor for /i = 0, that is when the front is straight and parallel to the the z-axis, 
h' = dh/dz is the spatial derivative of h and Ti. is the Hilbcrt transform of the derivative, 



(5) 
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Note that the perturbation expansion ([4]) of the stress intensity factor is invariant under translation of the front, 
/i — > /i + est. This expansion can be rewritten in the Fourier space as 

='5(9)-^l'?l^ + ^/ i2\qm-2f+eiq^e))h,h,^e^. (6) 



2. The equation of motion 

The equation of motion of the crack front can be obtained by considering the dynamic energy release rate for a crack 
propagating at velocity slow compared to the characteristic wave speed of the material. According to the generalized 
Griffith energy criterion [27i | , the crack must grow in such a way that the energy release rate G is equal to the dynamic 
fracture energy of the material T{v). This material property may depend on the instantaneous crack tip speed [23]; 
however, for low velocities, one can assume that the fracture energy is approximately constant T{v) ~ F. In the case 
of a propagating tensile crack front, this reduces to ^3] 



r, 



(7) 



where fi is the Lame shear coefficient and v is the instantaneous local normal velocity of the crack front. The velocity 
dependent function g{v), given in Appendix A, is universal in the sense that it does not depend on the details of the 
applied loading, nor on the configuration of the body being analyzed; it depends only on the properties of the material 
through the dilatational (cd) and the shear (cg) wave speeds. Last, K {t,v = 0) is the static stress intensity factor. A 
reformulation of Eq. ([7]) is the generalized Irwin criterion 28] , 



V^K [h] = , 

introducing the material toughness = \/2fiT. At low velocities, one can expand 

2 / / \ 3N 



hO - 

CR 'Z ' \CrJ \\CR 



(8) 



(9) 



where cr is the Rayleigh wave speed, a{K) is a dimensionless function of the material parameter k = (cd/cg) (see 
Appendix A). Eq. ([8]) now becomes 



V 

CR 



2 V CR 



1 K[h]-K, 
y{k) K [h] 



leading to the 2"'' order expression for the velocity 



V 

CR 



K[h]~K^ ^ 1 (K[h]-K^ 
" 2 



K[h] 



K[h] 



(10) 



(11) 



with a surprising simplification that makes the term inside the brackets independent of k. This implies the equation 
of motion 



1 



dh 



VTTh^ dt 



K[h]-K, 1 n<[h]-K, 
K[h] 2\ K[h] 



H {K [h] - 



[Fracture] , 



(12) 



where the Heaviside function H{-) accounts for the irreversibility of crack opening, the term vT+Ti^ is the ratio 
between the x-velocity and the normal velocity of the crack front, and U{k) = cr/q!(k) is a characteristic velocity. 
Note that C/(k), plotted in Fig. [21 is of the order of magnitude of the Rayleigh wave speed cr, which is very high. 
This gives a possible explanation for why even very close to the depinning threshold the movements of the interface 
are fast and therefore difficult to capture. 




K 

FIG. 2: The characteristic velocity U{n) in the equation of motion (|12|) . measured in units of the Rayleigh wave speed cr. 



B. Wetting 



1. The contact angle 



Consider a liquid film on a solid substrate, such that there is a triple contact line liquid/solid/gas where the thickness 
Y{x, z, t) of the liquid layer vanishes. The position x = /i(z, t) gives the location of the contact line (Fig. [3]). In order 
to establish its equation of motion, we first compute the contact angle as a function of h. We assume the motion 
of the contact line to be quasi-static. The length scales of interest are below the capillary length so that we neglect 
gravity. Moreover, the slope of the liquid layer is considered to be small. Within this framework, surface tension 
is solely involved in the shape of the layer and the thickness Y{x,z,t) is a solution of the Laplace equation. The 
experimental device imposes the slope of the film 9o far from the contact line. As a result, the problem is to solve 

( / \ 

< Y{x = h{z,t),z)^0, (13) 
dV 

— (x = -00,2;) = -6*0, 
V ax 

and more precisely to find the value of the contact angle 9 — —dV/dn at the location {x — h{z,t), z). 




FIG. 3; Schematic representation of a rough contact line. The straight reference line in the z-direction and the perturbation 
h{z, t) around it are also shown. 

This problem is the same as the one solved in [29l Isoj in the context of the diffusion-limited evaporation of a thin 
film. Translating into the terms of wetting contact line problem, this gives 

^ = 1 + 2H[h] -f \h''^ + hh" + mlh-HM - h'^n[h] + ih^nh)" + n ih^h") + mlh-Hlh-HM] + 0{h'^) . (14) 

At first and second order in h, the contact angle has the same terms as the stress intensity factor but with different 
numerical prefactors. This form was previously derived at second order in p7| . Note that the derivation of this 
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expression in [2^ [s^l suggests that higher order terms are polynomial in h''^ and in the operator T-i[h'-]. As will 
be discussed in the conclusion, this remark has implications on the relevance of high order terms in the equation of 
motion of the contact line. Eventually, for completeness we give the Fourier transform of the contact angle, 

+ / |-^|m|£(g - m - £) + ig^|m| + i|g|m^ - |g||TO||m + £|| . (15) 



2. Equations of motion 

According to the underlying microscopic physics, various equations of motion were proposed. In the model by 
Blake j3l| . the normal velocity derives from the surface energy of the system 



^ = uVl + h^ {0[h]' - 0l^) H(e[h] - , [Blake] (16) 



where U is the characteristic hopping velocity of the liquid molecules at the contact-line. The term VT + Ti^ is 
again the ratio between the x-velocity and the normal velocity of the line. 0eq is the equilibrium contact angle. The 
Heaviside H function might be used to account for the hysteresis in contact angle; 0cq then stands for the advancing 
contact angle. 

In de Gennes (3^ . the viscous dissipation is balanced by the the gain in surface energy 

rih 

— ^uVl + h^O[h]{e[hf-el^)H{e[h]~e,^). [de Gennes] (17) 

The capillary velocity U = 7//! is the ratio of the surface tension 7 o f the liquid and its viscosity ^. 
In the purely hydrodynamical approach due to Cox-Voinov [3^ . Is^ , 



— = uVl + N^ {e[hf - el^ H{9[h] ~ 6'oq) , [Cox-Voinov] (18) 

where U is again the capillary velocity. Note that the model with a diffuse interface proposed by Pomeau and 
Pismen (ssl ISo ] yields a similar macroscopic form. 



C. The annealed noise approximation 



In the case of propagation of a crack front through a heterogeneous material or wetting of a disordered substrate, 
the standard way to introduce the effect of heterogeneity in the equations of motion consists of modeling the material 
toughness Kc{z, h) or the equilibrium contact angle 9c^{z, h) as random functions which are separated into a constant 
and a fluctuating part as 

Kc{z,h) = K* + K*f]{z,h) , (19) 
6,^{z,h) = e* +6*f]{z,h) , (20) 

where r){z^ h) is a random noise with a zero mean and is assumed to be short-range correlated. 

Our approach is based on the assumption that one can approximate the noise term for the moving front, where 
h ~ Vt, by f]{x,h) ~ f](x,Vt) = r](x,t). This approximation amounts to saying that the fluctuations in the front 
velocity are small compared to the imposed driving velocity V. Also, we do keep nonlinear terms, since we claim 
(and will justify later) that they play an important role in roughening the interface. Obviously a linear equation of 
the kind described above, i.e. taking into consideration only the linear term in K[h] (or in 0[/i]), would not yield 
any roughness, and actually even if the KPZ nonlinearities (i.e. h'^ terms) are kept, we would also end up with a 
smooth surface, or at most logarithmically rough. This is a special case of the so called Fractal KPZ equation studied 
previously in [37| . 

The different equations of motion ()12|16|17|18|) . for fracture and wetting respectively, have the same mathematical 
structure, so that we can study both systems simultaneously. To this purpose, we rescale for fracture t — > t/U and for 
wetting h ^ h/A and t — > i/(4C/). When inserting the form K[h] (or 0[h]) in each equation, expanding to 2"'' order 
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in h and assuming the noise amplitude to be of the same order as h^, all the equations can be recast in the general 
form 



^ = w + n[h] - ^h'^ + ^hh" + H[hn[h]] + ^ {n[h]f + tj (z, t) , 

where the dimensionless velocity of the unperturbed front v and the coefficients a and b are given by 
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(21) 
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Resulting from Equations (|12I16I17I18P respectively. A comparison between these different cases is given below in 
Fig. m Notice that there is a real difference between the prefactors in the different cases in terms of sign, existence of 
zero-crossing, and general behavior as function of K* /Kq (or 9* /9o). 






K /K^ or 9 /e^ 



K /Kg or e /e^ 



K /K(j or 9 /Gp 



FIG. 4: The prefactors of the nonlinear terms in the equation of motion (|2ip as functions of the ratio of the average toughness 
to the driving stress intensity factor K* /Kq (or the ratio of the average equilibrium contact angle to the driving angle ^oq/^o 
for wetting). 



Interestingly, the velocity of the unperturbed front v in Eq. (|2ip can be eliminated, by transforming into a co-moving 
coordinate system, i.e. h h + vt, and so in the following we omit it. Finally, the noise correlations are described by 



{t] {z, t) rj (z', t')) = 2DS {z ~ z') 6 [t - t') , 



(26) 



where D is the variance of the noise. 



III. THE SELF CONSISTENT EXPANSION 
A. Description of the method 



The SCE method is based on going over from the Fourier transform of the equation in Langevin form to a Fokker- 
Planck form and on constructing a self-consistent expansion of the distribution of the field concerned. We thus 
consider the equation of motion (|2ip (transformed in the co-moving frame) in Fourier components 



dh. 



- -c h 

at ~ ' ' 



Mqimhihrn + Vq (0 



(27) 



where 



|g| and 



M, 



qim 



4v^ 



[\q\ 



a£m-b\£\ \m\ 



(28) 
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Q. being the linear length of the front. Note that in contrast to the KPZ problem Mqira has only the symmetry 
Mqim = M-q.-f..-m- Here, riq{t) is a noise term with zero average described by its variance 

{Vq (t) V,' it')) - 2DS {q + q') S it - t') . (29) 
Rewriting Eq. (|27p in a Fokker-Planck form, one gets 



dP_ ^ d 



dha 



D 



d 



jhq + ^ Mqirahlhr, 



P = 0, 



(30) 



where P({/iq},t) is the probability functional for having a height configuration {/ig} at time t. 

The expansion is formulated in terms of the steady-state structure factor (also called the two-point function) 

0, = {h^qhq) , (31) 

and its corresponding steady-state decay rate LOq, which describes the rate of decay of a disturbance of wave vector q 
in steady state and is defined by 



_,_ j;r{h-qmq{t))dt 

(h^.hq) 

From the scaling form ([2|), it follows that for small q^s, (j)q and LOq behave as power laws in namely 

u:q = B\q\\ 

where z is the dynamic exponent, and V is related to the roughness exponent by C = (F — l)/2 [ssj . 
The main idea of SCE is to write the Fokker-Planck equation dP/dt = OP in the form 



dP 



0(0) + + o(2) 



P. 



(32) 



(33) 
(34) 



(35) 



where Oq, Oi and O2 are zero, first and second order operators in some parameter. The evolution operator Oq is 
chosen to have a simple form 



d 



^qhq j , 



(36) 



whose corresponding zeroth order solution for P^^^ is a Gaussian given by 

pW({M)«exp „1^^ 



(37) 



implying that at leading order the two point function is given by = (f)q. Note however that (f)q and ujq 

arc still unknown functions. Next, an equation for the two-point function is obtained. The expansion has the form 

09 = 0(j + |0p,a)p|, where Cq is a functional of all 0's and cj's. This reflects the fact that the lowest order in the 
expansion is exactly the unknown (pq . In the same way, an expansion for ujq is given by cuq — LUq + dq 



|0p,a)p|. 



The 



two-point function and the characteristic frequency are determined by the requirement that the corrections to the 
zero order vanish, and thus given by the two coupled equations 



|0p,a)p| = 0, 
|0p,Wp| = 0. 



(38) 
(39) 



Eqs. (|38l39p can be solved exactly in the asymptotic limit of small g's yielding the required scaling exponents governing 
the steady-state behavior and the time evolution. 
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B. Equations for the structure factors 



Since the symmetries of the kernel Mg^mas defined in Eq. ([28)1 are different from those of the KPZ equation, we 
cannot use the SCE equations obtained in [22| and we need to derive them for the present problem. We first rewrite 
the Fokker-Planck equation ([50]) along the lines presented above as 



dP d 



+ uJphp 



d 



+ {cp - ujp)hp 



F = 0, 



(40) 



where 4>p and Cjp are unknown functions. The M term in Eq. is considered as the operator O^^^ and the differences 
{D — LLip4>p) and {cp — Cbp) are considered as O'^^ In a similar way we can expand the probability functional 

P{{hq}, t) = + + + . . . ^ (41) 

from which follows an expansion for an expectation value for an arbitrary observable .7^ to a prescribed order n 



in) 



(42) 



In order to calculate such averages we multiply Eq. (|40p by JFd/ig}) and integrate with respect to hq^s. After some 
integrations by parts we get 



dt 



dh^pdhp 



ujp < hp 



dhpl 



Y^Mp 



D — UJri 



dh-pdhp 



(n-2) 



dhpl 



(n-l) 



"dhpl 



(«-2)- 



0. 



(43) 



where we have added the superscript (n) to keep track of the order of the quantities with respect to the perturbative 
structure. At first we are interested in steady state properties and therefore drop any time dependence from the last 
equation. Later on, when deriving the equation for the characteristic decay rate the time dependence will be retrieved. 
Then, being interested in the 2-point function, we insert T — ^hqh_q into Eq. P5)) . To zeroth order we get 



0jq(f)q - 0jq{h^qhq)^°^ = , 

from which follows {h^qhq)^^^ — (j)q. Similarly, considering J- ~ ^hqh^q to first order gives 



{h^qhq)^^^ = <f>q- -j- Mqlmjh^hjnh- 



\(0) 



(44) 



(45) 



In order to calculate the correction we need to know the average of the 3-point function evaluated at zeroth order. 
To achieve this we insert !F = hLhMh-q (with L + M = q) into Eq. (|43|) and consider zeroth order contributions. 
This gives exactly zero, and so {h_qhq)^^'' = {h_qhqY^\ Thus, in order to get subleading corrections to the 2-point 
function wee need to go to second order 



UOq {h-qhqy"^' = UJq {h-qkqy^' + D - Cq {h-qhq^' - ^ Mqlm {h^hmh- 



(46) 



As above, in order to calculate {hihmh-q)'^^^ we insert J- — hLhuh-q (with L + M — q) into Eq. (j43p and consider 
first order contributions. This gives 

/U U h \(1) _ 2't'm{Mi^_rn,q + Me^q^_jn)+<t>eiMm-e,q + Mjn.q.-e) (Af-g,-£ -m + M-g.-m.-f) 
{nitlyyih^q) — ~(Pq — ;— 7 ;— 7 — — 7 ■ (41'] 



UJe + UJ,n + LOq 

Plugging this result into Eq. and using the symmetry of Mqgm gives 



OJl + LU„i + i^q 



~ /h h \(2) ~ /U U \(0) , n /h h I Mqim (Mqim + Mqmi) (I>e4'n 

UJq {h^qhqY ' = Uq {h^qhqY ' + D - Cq {h-qhqf ' + > — ^ ~ \ ~ 

Uli + LOra + 



UJe+UJm+ l^q 



(48) 
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The last step in the Self-Consistent Expansion is to impose self-consistency on this second order expression. This 
means choosing 0^ such that the two-point function to second order would be exactly equal to the zeroth order result. 

This is achieved by identifying {h-qhq)^"^^ = {h-qhq)^^^ = (pq and so the rest, denoted eq{(j)p,ujp} in Eq. ([38]). is set to 
zero : 

D - C 4> + (j> + Mqmi) {Mm-l,q + Mm,q-i) <pl ^ M^i'm {Mgi^ + Mgrnl) 4>e4'm _ ^ .^g^ 

' , use + LUrn + i^q „ OJl + LOm + iOq 

In a similar way we can derive an equation for tOq corresponding to Eq. ()39|) . Since tOq is related to dynamical properties 
(see Eq. we need to consider time-dependent quantities. Therefore, we can plug T = h^q{0)hq{t) into Eq. 

and proceed as above. As shown in [2^. [23j. an alternative way is to use Herring's consistency equation [30| which 
reads 

(Mqim + Mqml) {M,n.-l,q + Mm.q.-f) 0f _ „ ,f.„^ 

The last two equations allow us to determine the two unknown function (j)q and Uq. But before proceeding it is 
useful to rewrite them in the following way 

D - k|0g+/l(<Z)<^g+/2(<?) =0, (51) 

u;q - |<z| + J(q) = 0, (52) 



with 



^ J_ f [\q\(\e\+\q-e\)-2aiiq-e)+2b\e\\q-e\-e-'-(q-e)-']^^^^ 

^ J_ r [\q\{\e\ + \q-e\]-2a£(q-i)+2b\e\\q-i\-i''-(q-if][\q-me\MQn^^ ( rr^ 



327r J ^i+^g-i 

Finally, it is interesting to mention here that Eq. (|5ip can be understood as emanating from the short time balance 
of the original equation, while Eq. (|52p comes from its long time balance [1^. 



C. General solutions for the scaling exponents 

Eqs. (|51|52p can be solved exactly in the asymptotic limit (i.e. for small g's) to yield the required scaling exponents 
governing the steady-state behavior and the time evolution. The difficulty here arises from the fact that the integrals 
involved, Ii{q), ^2(9), and J{q)^ have contributions from large £'s as well as from small £'s. Therefore, one must 
consider the contribution of the large i integration domain on the small q behavior of the integrals (|53ll55p . For 
this, we break up the integrals Ii{q) and J{q) into the sum of two contributions If{q) + I^{q), and J'^{q) + J^{q), 
corresponding to domains of integration over f < A or £ > A respectively. We expand (q) and (q) for small q's 
and obtain the leading small q behavior of the integrals 

I>{q) = A,\q\~B,\q\cUq + Ciq\ (56) 

I>{q) = A2+B2 \q\ + C2q^ - D2 \q\ Uq , (57) 

J>{q) = A3 |g| - i?3 kl + C3g^ (58) 
where the leading order coefficients Ai,A2, A3 are 

(a + b~l){2b+l) f°° .(he 
Ai=A3 = ^ '- / e^^di, (59) 

Ibn Jj^ cue 

and so their exact values generally depend on the cutoff A. Also, note that ,63 = (it was kept for clarity). Using 
these results, Eqs. (I51l52p reduce to 

D + A2-{1- Ai) \q\ 4>q + /i< (<z) <i>q + (<z) = , (61) 
ujq - (1-Ai)|g| + J<(q) = 0. (62) 
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The advantage of Eqs. (|61l62p over Eqs. (|51l52p is that at the mere price of renormahzing some constants in both 
equations, we are left with the integrals Ii{q), ^^(9) and J^(g) that can be calculated explicitly for small g's since 
the power-law form for 0^ and oj^ for small £'s p3|34p can be used. The treatment of Eqs. (|61|62p is carried on by 
studying the various possibilities of balancing the dominant order for small q. Note also that the small g-dependence 
of each of the integrals (q) and J< (q) depends on the convergence of the integrals without cutoffs. So, to leading 
order in q one gets 

' \ F,\q\^-^-^ 4-r-z<o ' ^^^^ 

T<( \ / -^2 5-2r-z > 

^2 (9) - |F2|g|5-2r-. 5-2r-z<0 ' (^4) 
I<(a) ~ i ^il'^l 4-r-z>0 




FIG. 5: The possible domains for the exponents. Definition of the different sectors corresponding to the division of the (F, z)- 
plane by the curves 4 — F — 2 = and 5 — 2F — z = 0, with F > and z > 0. Note that the physical domain of the roughness 
of the front is given by 1 < F < 3, and thus solutions in sector 7 are not allowed. 

Let us consider the quadrant of the (F, z)-plane defined by F > and z > 0, where solutions may be expected. The 
lines 4 — F — z = and 5 — 2F — z = divide this quadrant into four sectors (see Fig. [5|) . The classical method [2^ 
is to investigate each sector separately and to decide whether or not a solution might exist there. In sector a defined 
by F < (5 - z)/2 and F < (4 - z), one can rewrite Eqs. (I61I62|) as 

D + A2 + E2-A{l-Ai-Ei)\q\^-'^ = 0, (66) 
B\q\^-il-A^-E^)\q\ = 0. (67) 

Eqs. (|66I67|) are self-consistent in the small q limit only if F = z = 1. These values are allowed by the defining 
conditions of this sector, thus they are possible solutions. In sector f] defined by (5 — z:)/2 < F < (4— z), the difference 
from sector a lies in Eq. (|66p which is now rewritten as 

D + A2-A{l-Ai~Ei) + F2\qf-^^-' = , (68) 

while Eq. ([67|) does not change and thus one still has z — 1. However, F = 1 is no more a possible solution of Eq. (|68p . 
since it is inconsistent with the defining conditions of this sector. The remaining option is given by 1 — F = 5 — 2F — z, 
which implies F = 3. This solution corresponds to a roughness exponent of C = 1 and thus is inconsistent with the 
assumption of small gradients used in the derivation of the equations of motion. Therefore, possible solutions do not 
exist in sector p. In sector 7 defined by (4 — z) < F < (5 — z)/2, Fig. [5] readily shows that one always has F < 1 and 
thus physical solutions can not be found in this sector. Finally, in sector 5 defined by F > (4 — z) and F > (5 — z)/2, 
Eqs. (|61l62p are rewritten as 

D + A2-A{l-A,)\q\'-^ + iAF,+F2)\q\'-'^-' = 0, (69) 
B\q\^-{1-A,)\q\+F,\q\'-^-^ = 0. (70) 

From the defining conditions of this sector the last term in Eq. (jTO]) is dominant over the second, which results in 
the scaling relation z = (5 — r)/2. Plugging this relation into the defining condition of this sector, one finds that 
one would necessarily have F > 3 (or equivalently C > 1), which is again inconsistent with the assumption of small 
gradients used to derive the equation of motion. 
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D. The rough phase 



Apparently, the only possibility one gets is the simple case T = z = 1 which corresponds to a moving flat front 
where perturbations propagate in a wave-like manner. However, a more careful inspection shows that one possibility 
was ignored, namely that of getting a fine-tuned case where the nonlocal elastic term is renormalized by the nonlinear 
term such that the pre-factor of l?]^"'" vanishes in Eqs (|66I68I69[) . This means that one has to re-perform the analysis 
by taking into account this possibility, for which the higher order corrections given in Eqs. (|56ll58p become important. 
In a previous publication ^20j , the classical analysis in the different sectors has been repeated while taking into account 
this new instability. Here, we adopt a more direct approach, namely reconsidering the original SCE equations ([5T|) - 
([52)1 . This approach clarifies the necessary conditions for the appearance of a rough phase. The idea is to expand 
Ii{q), hiq) and J{q) in powers of q for small q's 

oo 



Wtt J ujg 



, (a + 5-l)(2a^l) ^ iM^^^ ^ (l + 2a + 26-2«6-2.^) ^ , ^ 3) ^ (7,) 

CtlTX J LUf lOTT J UJg 



<^±t^/M., + OM), (72, 



16n J UJg 







Interestingly J{q) has the same expression as Ii{q) up to 2"'' order, which is the order that interests us here. When 
plugging this into Eqs. ([5T|) - ([5^ we get 



167r J LOg \ 167r J uog 

\ 



(a + 6 - l)(2a - 1) f e.',^g_^^ ^ (l-K2« + 26-2a6-2a^) ^ ^ ^^^^ 



1 




327r J ujf 167r 



00 

(a + 6- 1)(2&+ 1) [ l^<i)g 
16^? 




(a + 6-l)(2a-l) [ tM<t>t ^ {l + 2a + 2h - 2ab - 2a') ? ^^0. | _ ^ _ ^^^^ 



327r J ojj 167r J ujg 



12(0) only renormalizes the noise amplitude D and it is always a positive contribution. The analysis of the previous 
section is completely changed whenever the prefactor of |g| in the two equations above vanishes, i.e. whenever 

(a + b-mb+l) 7^^^..,,_i^o. (75) 



IGtt J UJg 



In this case the next order terms become important, and it immediately follows that oc q^' and ujq cx q'. Using 
previous notations, this means that F = z = 2 and Q = 1/2. Note however that in order to obtain this scenario a 
necessary condition is that the quantity (a -f & — 1)(26 + 1) must be positive. In Fig. [51 we plot this prefactor for 
the four cases discussed in this paper. Interestingly, in the case of wetting this prefactor is indeed positive for de 
Gennes' fEq. fTT)) and Cox-Voinov's (Eq. [T5)) mobility laws. However for Blake's mobility law (Eq. (|16p ). this prefactor 
is negative so that no rough phase can exist. Also, in the case of fracture the propagation law (Eq. [T^ . imposes 
that {a + h — 1)(26 -f l)/(167r) decreases with Kg/K* and even becomes negative for small Kq/K*, and hence with 
increasing mean velocity v of the front; the rough phase might disappear at high velocity of the elastic line. 

Actually, since the possibility F = z = 2 lies exactly on the boundary between the two sectors /3 and 5 (see Fig. [5|) , 
integrals such as in Eq. ([75]) diverge logarithmically for small q. So the scenario must be slightly corrected as follows 
(see appendix B for details). The branch described by F = z = 1 (the linear regime) rules for small q, q < qc- Then 
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FIG. 6: The existence of the rough phase is controlled by Eq. (|75|) and depends on the prefactor (a + & — 1) (2b +1)/ (IBtt) as a 
function of the dimensionless ratio K* / Ko (or 8cq/6o for wetting) for the four different cases discussed in this paper. As can 
be seen, in three out of the four cases, the positivity of this prefactor is satisfied for low velocities, where its value is around 
0.5. For Blake's mobility law, this prefactor is always negative (yet very small) and for fracture it becomes negative at high 
velocities. 



from Qc up to A (the upper cutoff introduced before), the scaling regime with F = 2 = 2 is dominant. The transition 
point is given by 

Qc ^ Ae"^ (0+6-1^x26+1) , (76) 

where A and B are defined by Eqs. (|33 p - (l34|) . The cutoff Qc is small as long as {a + b — 1)(26 + 1) > 0, which is the 
same condition as that encountered above, so that the rough behaviour exists on a wide range qc < q < A. This 
justifies the statement that the opposite signs in the prefactor of \q\ in the equations can drive the spectrum towards 
a rough regime. Moreover, it is now clear that the requirement of Eq. (j75p is not a strict one. It is sufRcient to make 
this difference very small in order to get a small qc (in other words gc ~* when the difference tends to zero) such 
that the rough branch of the spectrum takes over. We would expect then a structure factor of the form 

= I 2 ■ (77) 

qc\q\ + r 



IV. DISCUSSION 



In this paper we derived from first principles the nonlinear equations of motion for an in-plane crack front, and for 
a wetting contact line using the three available local velocity laws. These equations (I12I16I17I18P coupled with (|4ll4p 
could serve as reference for future research. 

Then we analyzed these equations with the aim of studying possible roughening of the front. To this end, we 
used the self-consistent expansion developed in [2^ and found the possibility of having a rough moving phase with 
a roughness exponent — 1/2 and a dynamic exponent z = 2. These results are relevant when Kq ^ K* (or 
9^9*) and are due to destabilization of the nonlocal elasticity by the nonlinear terms. Including higher order terms 
would not affect these results as they have the same structure as the quadratic nonlinear terms which we took into 
account (see Eg. Interestingly, since this destabilization is generic, it can be captured by other methods such as 
Dynamic Renormalization Group (DRG) used in p7j . For example, using our notations, the 1-loop expansion of the 
full propagator (corresponding to Eq. (61) in [13]) can be written as 

which resembles our result in Eq. ((73|) . having the same constant in front of the integral. Again, q^ relaxations in the 
system appear when the Igj-behaviour is killed by the nonlinear terms. 

We found that this scenario was consistent with three out of the four equations of motions we considered, namely the 
equation for the crack front and the two equations for the wetting front using the mobility laws of de-Gennes p7|) 
and Cox-Voinov (fT8|l . However, Blake's mobility law_(ll6|) was found to exclude the rough phase. This finding could 
serve as an argument against Blake's mobility law j3l| in the ongoing debate over the dynamical laws governing the 
spreading of fluids on solid substrates [3^ . 
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Our finding is consistent with the roughness exponents measured in crack propagation, ( ~ 0.6 [ll|, [iJl, in wa- 
ter/glycerol contact lines, C ~ 0.51 |15||. Our dynamic rough phase would apply to the quick motion of the lines in 
experiments. More precisely, a physical picture can be drawn as follows. In experiments, the driving force is not 
imposed at the front but mediated by a spring - the meniscus for contact lines and the elastic medium for crack 
fronts. As a consequence, the force acting on the front varies in space and time; it oscillates below and above the 
local threshold for motion. Above this threshold, the front moves quickly (an avalanche occurs) and its roughness 
corresponds to our calculation. Then the force drops below the threshold and the front freezes. Here comes into 
play the discontinuity in the equation of motion (the Heaviside function) corresponding to the irreversibility of crack 
opening and to the hysteresis in equilibrium contact angle. This discontinuity becomes important in the last steps 
of freezing and tends to further roughen the line, with a roughness exponent slightly larger than 1/2. Moreover, the 
roughness exponent should be larger for the very irreversible crack propagation than for the slightly irreversible wet- 
ting, as seen in experiments. To summarize, we propose that observations can be explained by a frozen dynamically 
rough interface. 

Interestingly, a similar phenomenon is implicitly present in the KPZ system. It is well known that any rough surface 
would eventually flatten by the KPZ system if the noise is stopped [Soj. However in real situations, it is compensated 
by a non-zero "angle of repose" that eventually freezes the system in a rough phase [41..] (this is expressed by an 
additional Heaviside function in the KPZ equation). It is also shown that the roughness exponent would be the same 
as that of the driven system if the freezing is done adiabatically [4l| . This shows that a hysteretic effect (existence 
of an angle of repose) can be consistently neglected once an above-threshold driving noise is present. It also hints 
that the final roughness of our system might be "history dependent" i.e. it might depend on the protocol of the 
loading/freezing, if not done adiabatically. 

Obviously, our study calls for more experimental and numerical work. In order to test our predictions, simulations 
such as in [21| could be performed for the nonlinear equations derived here and for a finite velocity of the elastic line. 
On the experimental side, several features could be tested: (i) the history dependance of the roughness exponent; (ii) 
the use of the form ([77| to fit structure factors; (iii) the transition from a rough phase (roughness exponent around 
1/2) to a flat phase at high velocity. Another experimental challenge would be the direct measurement of the dynamic 
exponent z predicted to be 2 here, which would provide an additional test. An interesting possibility would be to 
investigate the relaxation of the contact line after applying spatially periodic perturbations at given wavelengths as 
done in [i^]. 
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Appendix A 



In this appendix, we recall some elements of the elastodynamic theory of cracks [27| and apply them to the motion 
of slow cracks. We begin with the definition of the Rayleigh wave speed cr, which is the root of D{v) given by 



with as,d = Jl 



(79) 
(80) 

where Cd and Cg are the dilatational and shear wave speeds. In the expression for the dynamic energy release rate G 
in Eq. ([7|), we used the universal function [23] 



(l-as)arf(l-w/cR)^ 
D(v) l-v/cd 



2 /■ 
- — / arctan 

TT 



Ar]'^y/{ri'^ - 1) (k - 1]'^) 
(k-2?72)2 



(81) 



where k, = {cd/cg)'^. To derive the nonlinear equation of motion, we need an expansion of the hmction g{v) to 
second order in w, which can be written in full generality as 



,/^^l.a{K) — +piK)( — ) + 

CR V CR / 



O 



(82) 
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After algebraic manipulations, it can be shown that the coefficients a^n) and can be explicitly written as 



a (k) = 1 — —I arctan 

2 Crf Cd ' 



1 



dr] 

TT 



Pin) EE -a\K)., 



(83) 
(84) 



which in turn yields Eq. 



Appendix B 

In this appendix we find an estimation of , the transition scale between the rough branch described hy T = z = 2 
and the (logarithmically rough) linear branch described by F = z = 1. Our starting point are Eqs. (|73ll74|) . from which 
we get 



{\-W^)\q\+WM? 

|2 



(85) 



with 



u;, = (l-W^i)|(Z|+W^2|gr, (86) 



6.Z,+ (i±|^/M*. ,87, 



oo 

IGtt J uJi 



oo 

(a + fe- 1)(1 - 2a) f l'^Lo[(j}i (l + 2a + 26-2a6-2a2) f f-cf)/. 



W2 EE / -^dl~'— — '- / -^di. (89) 



From the equations above, we expect that at low g, Uq ~ \q\ and (pq — D\q\^^, which corresponds to the linear branch; 
whereas at larger q, ujq = -B|gp and (f>q — ^jg]"^, which corresponds to the rough branch (see Eqs. [33ll34|l . The 
transition between the two behaviors occurs when q = qc determined by 1 — Wi — W2qc- Using the small scale cutoff 
defined by A, which was introduced in Eqs. (|59l60p . we can estimate 

oo gc A oo 

'^d^= f'^di+f'-^di+f'^di (90) 
J uji J uji J LJe 





9c 




?c A 


OO 




de f 








-*j 


qc 


A 



^''^'d£^^Dq'^ + ^log-+ I ^-^d£, (91) 



ujd ^ B qc J uji 

A 

and deal with similar integrals in the same way. Hence the defining equation of qc becomes 

(a + 6- 1) (2a- 1) 




(92) 
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In the limit Wi — > 0, we get — > 0, so that we can keep the most dominant terms and get 



qc= A exp 



( 



B IGtt (a + 6 + 2) 



A{a + b-l){2b+l) (a + 6- 1)(26+ 1) 



) 



(93) 
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